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^ Abstract. In this paper, we prove that in small parameter regions, 

arbitrary unitary matrix integrals converge in the large limit and 
match their formal expansion. Secondly we give a combinatorial model 
, for our matrix integral asymptotics and investigate examples related to 

free probability and the HCIZ integral. Our convergence result also leads 
us to new results of smoothness of microstates. We finally generalize our 
, approach to integrals over the othogonal group. 
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Introduction 



Matrix integrals provide models for physical systems (2D quantum grav- 
ed ' itation, gauge theory, renormalization, etc.), and generating series for a 
^ ■ wide family of combinatorial objects (see e.g pOl [28] ) . 

Q>^ I Gaussian integrals are the most studied. It was shown by Brezin, Itzyk- 

son, Parisi and Zuber |7| that perturbations of Gaussian integrals expand 
, formally as a generating function of maps, sorted by their genus when the 

I dimension of the matrices is regarded as a parameter. Such 'topological' 

■ expansions were shown also to hold in the large N limit, and then to match 

! with the formal expansion on a mathematical level of rigor by two authors 

I [m [iTl [23] and previously in the one matrix case in [H [2] and [12] . The 

relation of Gaussian matrices with the enumeration of maps is an easy con- 
sequence of Wick calculus -or equivalently Feynman diagrams- see [28] for 
a good introduction. According to 't Hooft [20j, such topological expansion 
^ , should hold in the more general context of models invariant under unitary 

5-H I conjugation. This leads us to concentrate in this article on matrix integrals 

^ given by 

(1) In{V, Af) := [ e^Tr(y(c/„c/;,Af ,i<i<m))^f^^ . . . 

where {A^ , 1 < ^ < rn) are N x N deterministic uniformly bounded matri- 
ces, dU denotes the Haar measure on the unitary group (normalized so 
that J^^ dU = 1) and V is a polynomial function in the non-commutative 

variables {Ui,U* , Af ,1 < i < m). Tr denotes the usual trace on x 
matrices given by Tr{A) = Yli=i ^a- 

We will study in this article the first order asymptotics of matrix integrals 
given by ^ when the joint distribution of the (Af , 1 < i < m) converges; 
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namely for all polynomial function P mm non-commutative indeterminates 

(2) lim ^TriP{Af, 1 < i < m)) = t{P) 

for some linear functional r on the set of polynomials. Without loss of 
generality, we will assume that {Af, l<i< m) are Hermitian matrices. 

For technical reasons, it is convenient to assume that the polynomial V is 
such that Tv{V{Ui, U*,Af , 1 < i < m)) is real for all Ui G Un, all Hermitian 
matrices Af , for all i G {1, • • • , m} and G N . 

Under those very general assumptions, the only result proved so far is the 
formal convergence of these matrix integrals. Namely, it was proved in [8] 
by one author that for each k, the quantity 

l^iV-^ log / e^^^<''^^-^tA-^<^<^))aU, . . . dUrnl=, 

converges towards a constant fk{V, r) depending only on the limiting distri- 
bution of the Af's and V. Besides, if V is polynomial with integer coeffi- 
cients, fk{y, t) is a polynomial function with integer coefficients of the limit 
moments of the A^^s. 

In this paper we will answer affirmatively to the following, previously 
open questions: 

(1) Does the limit of the matrix integrals exist for small parameters z? 

(2) Does the power series Ylk Ij/fcl^)''") have a strictly positive radius 
of convergence? 

(3) Is the limit of the matrix integral equal to the power series? 
The following Theorem is a precise decription of our results: 

Theorem 0.1. Under the above hypotheses and if we further assume that 
the spectral radius of the matrices {Af ,\ < i < m,N ^ N) is uniformly 
bounded (by say M), there exists e = e{M,V) > so that for z G [— e,e], 
the limit 

FvAz) := hm 4 l°g / e^^TriViU,,u:,Af^<^<m))^u^ , . , ^jj^ 

exists. Moreover, Fy^riz) is an analytic function of z £ Cn B{0,e) = {z £ 
C : \z\ < e} and for all /c G N, 



= fk{V,T). 

2=0 



This also implies that the series Fy^.,-{z) has a positive radius of conver- 
gence, a result which had not been proved by the techniques of [8j based on 
Weingarten functions. 

Our approach is based on non-commutative differential calculus (in par- 
ticular on the resulting Schwinger-Dyson or Master loop equations) and 
perturbation analysis as developed in the context of Gaussian matrices in 
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[Ml [m [23 j . Another possibility to prove the equaUty between real and for- 
mal limits would have been to show convergence of the integrals for complex 
parameters z. We have not yet been able to follow this line successfully, and 
this remains an open question. 

An important example of unitary matrix integral is the so-called spherical 
integral, studied by Harish-Chandra and by Itzykson and Zuber, 

HCIZ{A,B):= [ e^Tr(f/*At/B)^f^_ 
JugUn 

This integral is of fundamental importance in analytic Lie theory and was 
computed for the first time by Harish-Chandra in [T9]. In the last two 
decades it has also become an issue to study its large dimension asymptotics 

[m [Ml El US]. 

Theorem 10.11 holds true for the HCIZ integral. It thus relates the results 
of [8] (which computed the formal limit of the HCIZ integral) and those of 
|18j (where the limit of HCIZ{A, B) was obtained (regardless of any small 
parameters assumptions) by using large deviations techniques) . Let us recall 
the limit found in |18J. Let us define 

= i^lJ'ix^) + \ j y log |x - y\dn{x)dn{y). 

If (resp. hb) denote the limiting spectral measure of A (resp. B), assume 
that I{^jla) and lifJ-s) are finite. Then, the limit of N~'^ log HCIZ (A, B) is 
given, according to [12], by 

(3) 

/(^A,/Xij) := lim log HC I Z {A, B) 

Af— >oo iV^ 

where the infimum is taken over m, p so that the measure- valued process 
Utidx) = pt{x)dx is a continuous process, = pA, pi = PB and 

dtpt{x) + d^rritix) = 0. 

The inf over {pt,mt) is taken (see }14j) at the solution of an Euler equation 
for isentropic flow with negative pressure —^p^. 

Theorem ID. II shows that I{p p depends analytically on (5 in a 
real neightbourhood of 0, a result which is not obvious from formula ([3]). 
Moreover, the coefficients of this expansion count certain planar graphs (see 
section [5]), as summarized in the following theorem. 

Theorem 0.2. Denote \/^f)P the probability measure 

y/PHf) = I f{Vf^x)dp{x). 
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Assume that o-nd are two compactly supported probability measures. 
Then, there exists Po > such that for all (3 € [0, (3q\, 

n>0 

converges absolutely. Moreover, we have 

MnifXA, I^b) = ^ M^{ha,1^b)- 

m admissible maps of s„ 

Tin is the set of planar maps drawn above n vertices defined as stars of 
type U*AUB by gluing pairwise oriented arrows and possibly rings and 
MminA, IJ-b) is the weight of the map. 

We refer the reader to section [5] for the definitions of stars, admissible 
maps and weights. Our definition of planar maps is more complicated than 
those arising in the topological expansion of Gaussian matrix models (and 
which are directly related with Wick Gaussian calculus and Feynman dia- 
grams): the sums are signed and we have a notion of admissibility. However 
it was an open question in mathematical physics to have a graphical model 
for unitary integrals (see [36]). Moreover, this graphical interpretation gives 
a new understanding of cumulants formulae (see section 16. 2p . 

The convergence of other integrals was still unknown and it is one of 
the points of this paper to show their convergence. We use it to study 
Voiculescu's microstates entropy evaluated at a set of laws which are small 
perturbations of the law of free variables, and prove regularity of microstates 

Theorem 0.3. For tracial states n satisfying suitable assumptions described 
in Theorem \8.1\ 

X{^i) := liminfliminf -^log/u®" {TR{fx,e,k)) 
= lim sup hm sup log (T^ (// , e , /c) ) 

fc| oo 

and a formula for x(/^) can be given. 
This result generalizes section 4 in |16] . 

The paper is organized as follows: after setting our working framework 
(section[l]), we study the action of perturbations upon the integral lN{y, ^f) 
and deduce some properties of the related Gibbs measure; namely that the 
so-called empirical distribution of the matrices under this Gibbs measure 
satisfies asymptotically an equation called the Schwinger-Dyson equation 
(section [2]) . Then, we study this equation and obtain uniqueness for param- 
eters of the potential V small enough (section [3]) and analyticity (section 

We also describe a (new) combinatorial solution of Schwinger-Dyson equa- 
tion (section [5]) and therefore of the first order of unitary matrix integrals. 
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We deduce applications of these results to free probability (section [6]) and to 
the convergence of matrix integrals In{V, ^j^) (section [7]) . Finally, we point 
out some consequence of our result for free entropy (section [8]) . At last, we 
consider the case where integration holds over the orthogonal group instead 
of the unitary group and show that the first order of such integrals is the 
same. 



Let Un be the set oi N x N unitary matrices, A4n the set of N x N 
matrices with complex entries, TCn the subset of Hermitian matrices of A^at 
and An the subset of antihermitian matrices of A^at- We let m be a fixed 
integer number throughout this article. We denote by (^f^)i<i<m a m-tuple 
of NxN Hermitian matrices. We shall assume that the sequence (^i^)i<i<m 
is uniformly bounded for the operator norm, and without loss of generality 
that they are bounded by one. 



1.1. Free *-algebra. Let C((C/j, C/*, Aj)i<i<m) denotes the set of polyno- 
mial functions in the non-commutative indeterminates (t^2, ^i)i<2<m 
with the relation 



Note that in general we may want to consider models with a number of 
"deterministic" indeterminates Ai different from the number of "random 
unitary" indeterminates Ui, but this general case can be obtained from the 
previous one by looking only at a sub-algebra and our convention shortens 
a little the notations. The algebra C{{Ui,U* , Ai)i<i<m) is equipped with 
the involution * so that A* = Ai, (Ui)* = U*; (U*)* = Ui and for any 
Xi,--- ,Xn € {Ui, U* , Ai)i<i<rn, any z £ C, 



{2X1X2 ■ ■ ■ Xn^iXn)* — zX^X^^i ■ ■ ■ X2XI. 

Note that for any Ui £ Un, Ai £ TCn, and P G C{{Ui, U* , Ai)i<i<m) , 
{P{U, U:,A.„l<i< m))* = P*{Ui, U*,A^,l<i<m) 



where in the left hand side * denotes the standard involution onM^. We de- 
note C((;7i, C/*, Ai) 

i<i<m)sa the set of self-adjoint polynomials; P — P* , and 

C{{Ui,U*,Ai) 

i<i<m)a the set of anti-self-adjoint polynomials ; P* — —P. In 
the sequel, except when something different is explicitly assumed, we shall 
make the hypothesis that the potential V belongs to C{{Ui, U* , Ai)i<i<m)sa, 
which insures that Tr {V{{Ui,Ul , A'^)i<i<m)) is real-valued for all Ui G Un 
and E Hn- Conversely, any potential V such that Tr {V{{Ui, U*,Af))) 
is real-valued for all Ui G Un and Af G Hn is self-adjoint up to the addi- 
tion of some commutators (which does not change the trace). Indeed, this 
implies that Ti {{V - V*){{Ui,U* , Af^) 

i<i<rri)) vanishes, which insures that 
V — V* = '^[PiQi — QiPi for some polynomials Pi,Qi, cf [9] Lemma 2.9 



1. Notations 




UiU* = U*Ui = 1. 
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for a probabilistic proof or [22j, Proposition 2.3 for a direct proof (in the 
real symmetric case, but directly adaptable to the Hermitian case). Then, 
W:=V + EiiQiPi - PiQi)/'^ is self-adjoint. 

1.2. Non-commutative derivatives. On C{(Ui,U* , Ai)i<:i<m) , we define 
the non-commutative derivatives di, 1 < i < m, given by the linear form 
such that 

diAj = 0, diUj = U=jUj (S) I diU* = -li=jl(S)U*, Vj, 
and satisfying the Leibnitz rule, that for P,Q £ C((C/j, ?7*, ^i)i<j<m), 

(4) diiPQ) = diPx {l®Q) + iP^l)xdiQ. 

Here, x denotes the product Pi Qi x P2 <Ei Q2 = P1P2 ® QiQ2- We also 
let Di be the corresponding cyclic derivatives such that if m{A i?) = BA, 
Di = mo di. 

If g is a monomial in C{{Ui, U* , Ai)i<:i<m) , we more specifically have 

(5) dig = ^ qiUi®q2- ^ qi®U*q2 

q=qiUiq2 q=qiU*q2 

(6) Dig = ^ q2qiUi - ^ f//g2gi- 

q=qiUiq2 q=qiU*q2 

1.3. Bounded tracial states. Let T be the set of tracial states on the 
algebra generated by the variables {Ui,U* , Ai)i<i<m, i-e. the set of linear 
forms on C{{Ui, U* , Ai)i<i<rn) such that for all P,Q £ C{{Ui, U* , Ai)i<i<rn) , 

fi{PP*) > 0, fi{PQ) = f,{QP), /x(l) = 1. 

Throughout this article, we restrict ourselves to tracial states fi £ T such 
that 

;u((4*Ai)") < 1 Vn G N, Vi G {!,••• ,m}. 
We denote this subset of T. 

Note that for any monomial q G C((C/j, [/*, ^j)i<i<m), the Holder's in- 
equality implies that for any ^ £ 

(7) /i(gg*) < 1. 

We endow ^A with its weak topology: /x„ converges to // if and only if for 

ah P G C{{Ui,U*,Ai) 

l<i<m/i 

lim /i,(P) = ^(P). 

By equation ([7]) and since the above topology is the product topology, M. 
is a compact metric space by Banach Alaoglu's theorem. 

We denote fi^ the empirical distribution of matrices Af G Hn and Ui G 
Un which is given for all P G C{{Ui,U* , Ai)i<i<m) by 

= ^Tv {P{U„ U*,Af,l<i< m)) . 
This object will be of crucial interest for us. 
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The notation stands for the set of tracial states of M re- 

stricted to the algebra generated by the (j4j)i<j<m- In particular, the limit- 
ing distribution r given by ([2]) belongs to -^|(Ai)i<i<™- 

1.4. Tracial power states. Let V G C((f7j, C/*, Ai)i<i<m)sa and Hy be the 
distribution on given by 

fi^idUi, • • • , dU^) = In{V, exp(7VTr(y))dC/i • • • dUm- 

We define, for ah P £ C{{U^, U* , A,)i<,<^) , 

. r f ^TrPe^^^dUi-.-dUn 

' '^^^^ Je^^''ydUi...dUn 

In the following, an n-tuple of monomials {qi)i<i<n in C((C/j, U* , Ai)i<i<m) 
will be fixed and we shall take V = Vt = X^"=i Then, Jly^{P) can be 
seen as a power series in the ti's; 



n|k| 

(8) -^l{p)--=T. 



k! n, dt\^ 



We will call /U a 'tracial power state' of M if and only if it is a map 

l<i<m/ ^IIT'll 

with for all a, 6, fi{ab) = jji{ba). Here C[[t]] is the algebra of power series in 
the variables ti, • • • In particular, we may view as a tracial power 
state of M.. 

1.5. Cumulants. The classical cumulants {Cfc}fc>o are defined via their 
formal generating function: 



\ogE{e''')=Y,t''Ck{X,...,X)/k\ 



k>0 



This equality holds also in a complex neighborhood of for t if X is bounded. 
We also define the cumulants for k in N": 



logS(e*^^i+-+*-^") = t^Ck(Xi, . . . ,X„)/k 



I 

keN" 

k _ rr 



where k = (/ci, • • • , k! = H- ki\, |k| = ki and t"* = Hi V- Note that: 

Ckl-'^l, • • • , Xk) = C|k|(-^l, • • • ,Xi, - ■ ■ , Xn, • • • , -'^n) 

where in the previous list the variable Xi appears ki times. 
Let us recall some properties of these cumulants. 

Proposition 1.1. The following two statements hold true: 
(1) 

—^-^^——^-^= ^ t Ci,k(y,Xi,...,X„)/k! 

^ ^ keN" 
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(2) 

E(^QtlXl + ...+tnX„^ E(^^tlXl + ...+t„Xn'^ £](^f^tlXl + ...+t„X„-^ 

= 5^t'^Ci,i,k(y,^,^i,...,^n)/k! 

fc>0 

Proof. 

Item (1) is obtained by replacing tiXi+. . by yY+tiXi+. . 

and differentiating the generating function of the cumulants in y at y = 0. 

Item (2) is obtained by replacing tX by yY + zZ + tX and differentiating 
the equality defining the cumulants in y and z at y, z = 0. 

□ 

2. Matrix models 

We first investigate the asymptotic behavior of the random state /i^ under 
as a random tracial state. We then consider Hy = ^y{(j,^) evaluated 
at a polynomial and study its convergence as a formal power series in the 
parameters of the potential V . We show that they satisfy asymptotically the 
same type of equations called Schwinger-Dyson (or Master loop) equations. 

2.1. Behavior of (i^ . The main result of this section is the following 

Theorem 2.1. Assume that V is self-adjoint. For all polynomial 

l<i<m/ ! 

lim {fL^(g)fi^{diP) + fi^{DiVP)}=0 n^a.s. 

N—*oo 

In particular, any limit point fi & M. of fl^ under Hy satisfies the Schwinger- 
Dyson equation 

(9) ti ix{diP) + n{D,VP) = 

for all P G C{{Ui,U*,Ai)i<i<m) and ^|(a,)i<,<,„ = 

The idea of the proof, rather common in quantum field theory and success- 
fully used in |16 p i7 1 [23]. is to obtain equations on fi^ by performing an infini- 
tesimal change of variables in lN{y, ^f)- More precisely we make the change 
of variables V = (Ui,--- ,U^) £ ^ *(U) = (^'i(U),--- ,^'™(U)) G 
Wf} with 

where the Pj are antisymmetric polynomials (i.e. P* = —Pj)- This change 
of variables becomes very close to the identity as N goes to infinity, reason 
why it is called "infinitesimal" . 

Lemma 2.1. The function ^ is a local diffeomorphism and its Jacobian 
has the following expansion when N goes to infinity 

I det J^(U)| = e^^' rr®Tr(a,P,(C/.,c/*,Ai,i<i<m))+o(i) 
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where 0(1) is uniform on the unitary group (but may depend on P). 
Proof. 

Let us first recall the following two elementary results of differential ge- 
ometry: 

(1) The map exp : M.n — TWat is differentiable and: 

DiffM exp .H := lim £-i(e*^+^^ - e^) = T^^^l 



where AdM is the operator defined by AdMH = MH — HM. 

(2) If P E C{{Ui,U* , Ai)i<:i<:m) is considered as a function of the Ui's, 
then it is differentiable and its differential with respect to the i-th 
variable in the direction A, for A in An, is 

Di&iP.A := lime-^PiUi, •••,[/,-!, Uie'^, Ui+i,---)- P(U)) = diP^- 

£— >0 

As a consequence, if we fix A in and i € {!,••• , m}, one has 
Diffi*,(U).^ = li=jUjA + ?7,Diff A p.(f;) exp .(A^.p^.y^) 

X +°° (Ad A p. 

= li.jUjA + Uj^^,i(V)A. 
with $jj(U) the linear map from An into M.n given by 

+ 00 CA(J ^ 

We can factorize the term C/j to obtain 

(10) Diff^'(U) = U o {May^ + A$(U)) 

with U o (Ml, ■ ■ ■ , Mjn) = {UiMi, ■ ■ ■ , UmM^n) and $ the linear operator 
from A^ to A^jv whose blocks are the $ij(U). 

Since the operator norms of the Ais and the UiS are uniformly bounded 
in N , the operator norm of AdAp.(u) ^ operator on {M.n, \\-\\oo) is also 
bounded. Thus, $jj(U) is a uniformly bounded operator from An to AIat- 
Thus, the norm of -^^>(U) is less than 1/2 for N large enough. For those 
A^, ^' is a local diffeomorphism with positive eigenvalues. 

We can now compute the Jacobian of ^ 

|detJ*(U)| := |detDiff^'(U)| = | det C/|| det(/ + A$(u))|. 
It can be easily checked that | det U\ = l. 
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Besides, the positivity of the eigenvalues of / + A$(U)/iV allows us to 
replace the determinant by the exponential of a trace: 

I det MU)\ = exp(TVlog(7 + ^$(U))) = exp ( - J] il^TV($(U)f) 



P> 



^ pNP 



Note that since <1> is a bounded operator on An, which is a space of dimension 
N'^, the p-th term in the previous sum is at most of order N'^^'^\ We only 
look at the terms up to the order 0{N). A quick computation shows that if 

An An 

is considered as a real endomorphism, Tr(p = ^^Tr^^Tri?/. Indeed, if we 
consider E{kl), 1 < k,l < N the canonical basis of An, 

Tp(l.l\ . /-^ ^r=k,j=l + lr=l,j=k 



for < / and 



for k>l,TT^p = Tt{Eikl)*ip{E{kl))) = TiAiTtBi. This is sufficient 
to obtain the first term of the Jacobian: 



^IV($(U)) = -^Tr($,,(U)) = - J2^®TT{d,P,{Uj,U*,A,)) + 0{l) 

i i 

with 0(1) is uniformly bounded on U^. Here we used that the operator 
norm of Ad_x p.^^) uniformly small. 

□ 

Before making the change of variables we show that * is a bijection. 
Lemma 2.2. For N large enough, is a diffeomorphism ofUjlJ. 
Proof. 

First observe that since is a local diffeomorphism, its image is open 
in Z//™. Besides, since is compact and ^ is continuous, the image is 
compact and therefore closed. Thus by connectedness of UJ!}, and since 
'J'(ZY^) is closed, open and non-empty, ^' is surjective. 

The only property we still need to prove is the injectivity of ^. If "^{U) = 
^{V) then for ah j e {1, • • • , m}, 
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Thus, if is sufficiently large so that ■^Pj{U) is in a domain where the 
function exp is 2-Lipschitz, we obtain 

lie/, - F.lloo = \\U,V; - 1|U = ||e^^^(^)e^^^(^) - 1|U 

= ||e^^^(^) -e^^^^lU < ^Wm) - m)\\oo 

with 1 1. 1 loo the operator norm. Since (Pj, 1 < j < m) are uniformly Lipschitz 
on U^, we conclude that X^jLi II ~ ^ lloo vanishes for sufficiently large N. 

□ 

We can now prove Theorem 12. II 
Proof. 

Let us define 

^""(P) = E (l^MD^VP,) + ^Tr ^Tr(a,P,)) • 

i ' 

We expand Try(^'(U)j, ^'(U)*, Aj, 1 < i < m) as 

Tr(y(^'(U)„ ^(U)*, .4,, 1 < i < m)) - Tr(y(C/i, [/*, 1 < i < m)) 

and perform the change of variables U — > ^'(U) in Ijqiy^A^y, 

^ I g7vTr(y(*(U).,*(u)*,A.,i<.<™))| det J^(U)|dC/i • • • df/^ 

= y ^NY^(P)+Oil)^NTriViU.,U*,A.^<^<m))^JJ^ . . . ^jj^ 



where we used (llip and Lemma [2. 11 0(1) is of order one independently of N 
and uniformly on the unitary matrices (C/i, • • • , Um)- Thus we have proved 
that 

j e^^"(^)d^{y(U) = 0(1). 

Borel-Cantelli's lemma thus insures that 

limsupy^(P) < a.s. 

N^co 

and the converse inequality holds by changing P into —P since Y'^ is linear 
in P. This proves the first statement of Theorem 12.11 The last result is 
simply based on the compactness of M and the fact that any limit point 
must then satisfy the same asymptotic equations than jl^ . 

□ 
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Another consequence of this convergence is the existence of solutions to ([9]) 
for any self-adjoint potential V (since any limit point of fi^ in the compact 
metric space A4 will satisfy it) a fact already proved in [6]. Moreover, since 
these solutions are limit points of jl^ , they belong to Ai and in particular 
< 1 for any monomial q. 

2.2. Moments of fi^ . In the sequel, we denote by E the expectation with 
respect to the Haar measure on the unitary group. The goal of this section 
is to show (see Proposition [2TT]) that cumulants also satisfy a formal version 
of Schwinger-Dyson equation. We start with the following lemma: 

Lemma 2.3. One has, for all i all N, all monomials qi,--- ,qn and all 
k = (A;i,-- - ,kn) in N", 

N'E (a^ {d,P) . (A^(gi))'^ • • • (A^(g.))'") 

+ Y: k,E {{(."^{q,))'^ ■ ■ ■ if.'^iq,))'^-' ■ ■ ■ (A^(g„))'" A^(A<Z, ■P))=0 

j 

Proof. 

Following Lemma 12.11 we write down the change of variable 

: U ^ (C/i, • • • , C/i_i, ?7ie^^-(U), C/i+i, • ■■ ,Um) 

in the integral / {{fi^qi)^^ ■ ■ ■ {fi^ qn)^")dUi ■ ■ ■ dUm, where the integration 
is taken over the unitary Haar measure. Its Jacobian satisfies 

I det J^{U)\ = l + ^Tr ^ Tr{diP) + o(A). 
and we have the expansion 

TT{qj{^{V)i,^{V)*,Ai,l< i<m))= Tv{qj{Ui,U* , A,! < i < m)) 

+ \Tr{Diqj ■ P{Ui,U*,Ai, 1 < i < m)) + X^oiX) 

where the o(A)'s are for a given P uniform bounds in N. The first order of 
the Taylor expansion of this change of variables around A = proves the 
claim. 

□ 

Proposition 2.1. As a formal series equality, one has, for all i, for all t, 
E[f,^ ^ A^(9iP)e^''^"(^*)] + Elf.'^iD.Vt ■ P)e^'^"(^*)] = 0. 
Proof. 

Multiplying the equality of Lemma [2. 31 bv t''A^^I''l~^/k! and summing over 
k in N" gives the desired identity. 

□ 

Finally we study the large N limit /x-^ of these formal states (the index / 
stands for "formal"). 
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Theorem 2.2. Let Vt he the polynomial X]j=i^i9i- -^'^'^ -f; sequence 
jly^{P) converges as a formal series (i.e. coefficientwise) when N goes to 
infinity to some fi^{P). Besides, fj,-^ satisfies the family of equations, for all 
i, for all P, 

=0. 

Proof. 

First, we prove the existence of a hmit. By the first item of Proposition 
ll-H we can express fly^{P) as a sum over cumulants, 

Ayt(^) = E t'^Ci,k(^TrP,iVTr(?i,--- ,iVTrg„)/k!. 

keN" 

The Hmit in N, of the Ci_k(^TrP, A^Trgi, • • • , NTiq n) was proved to exists 
in [8] so that fi^ is well defined. 

Item (2) from Proposition 11.11 implies 

ij;(^TrPi^TrP2e^^^^) j^(^TrPie^Try) j^^^^^^^^Ttv^ 

E{e^TTV^ ^(e^Try) ^(e^TrV) 

= E n^i'i'l'^^'^'^i' ^TrP2,iVTVgi, • • • , iVTrg„). 

A:>0 

Now, it follows from [8] that elements on the right hand side have decay 
A'^"^ so that the coefficientwise limit is zero. This can be interpreted as a 
formal convergence of measure result for the states fi^ . 

The proof of the Theorem follows from this observation and from Propo- 
sition 12.11 

□ 

3. Study of the Schwinger-Dyson equation 

We have shown that the limit points of the matrix model satisfy the 
Schwinger-Dyson equation Q. The aim of this section is to study this 
equation and show that it has a unique solution. 

Definition 3.1. Let r G M.\( a.\,^.^ . A tracial state u £ A4 is said 
to satisfy Schwinger-Dyson equation SD[V,r] if and only if for all P E 

l<i<m/ ; 

KP) = r{P) 

and for all P £ C{{Ui, U* , Aj)i<j<m), all i £ {1, ■ ■ ■ , m}, 

fi ^ f^idiP) + n{D,V P) = 0. 
Let V £ U*,Ai) 

i<i<m) ■ One can decompose 1/ in a sum 

n 

V = Y,t,qi{Uj,U*,Aj,l<j<m) 

i=l 
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with monomial functions qi and complex numbers tj. The monomials {qi,l < 
i < n) will be fixed hereafter. We let D be the maximal degree of the 
monomials Qi. 

Here we prove that r is uniquely defined provided that the parameters 
(ti,l < i < m) are small enough. 

Theorem 3.1. Let D an integer and r a tracial state in -^|(yii)i<i<„ be 
given. There exists e = e{D,m) > such that if \ti\ < e, there exists at 
most one solution fi to SD[V,r]. 

From this and Theorem 12.11 we deduce the following 

Corollary 3.1. Assume that V is self-adjoint. Let D an integer and r a 
tracial state in -^|(^i)i<i<„ be given. There exists e = e{D,m) > such 
that if \ti\ < e, fi^ converges almost surely to the unique solution ix of the 
Schwinger- Dyson equation. Moreover, fty = ^y{jl^) converges as well to 
this solution as N goes to infinity. 

This result is obvious since Theorems 12.11 and 13.11 show that jl^ has a 
unique limit point, and thus converges almost surely. The convergence of Jly 
is then a direct consequence of bounded convergence theorem since fi^ E M. 

Actually Theorem 12.11 and Corollarv 13.11 do not use the assumption that 
the matrices {Af , 1 < i < m) are deterministic, but only that they are 
bounded and have a converging joint distribution. Therefore these two re- 
sults extend to the case where these matrices are random, independent of 
the {Ui,l < i < m), and satisfy the above two conditions almost surely. 
This observation implies that our result can also encompass the case of the 
truncated GUE or other classical bounded matrix models. 

We prove now Theorem 13.11 

Proof. 

Let /u be a solution to SD[V,t]. Note that if we take q a monomial in 
C{{Ui, U* , Aj)i<j<m), either q does not depend on {Uj, U* ,1 < j < m) and 
then ^{q) = T{q) is uniquely defined, or q can be written as q = qiUfq2 
for some i £ {I,-- - ,m}, a G {— 1,+1} and monomials qi,q2 (Here, note 
that = U*). Then, by the traciality assumption, fi{q) = n{q2qiUf ) = 
^{Ufq') with q' = q2qi- Remark that we can assume without loss of gener- 
ality that the last letter of q' is not Uf"". We next use SD[V,r] to compute 
li{Ufq) for some monomial q. We assume first that a = — 1. Then, by dH), 



^^ {U*q) = 



-l^{U*q) + U* ^Ixdiq. 
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Taking the expectation, we thus find by ([5|), since ^(1) = 1, that 

fi{U*q) = fi(g> fj,{U* (g> 1 X diq) + fi{DiVq) 

= KQi)f^{Q2) - fi{U*qi)fi{U*q2) 

q=qiUiq2 q=qiU*q2 

(12) 

j 

where {tij,qij) are such that DiV = ^jtijqij. Note that the sum runs at 
most on Dn terms and that all the tij are bounded by max|ij|. A similar 
formula is found when a = +1 by differentiating qUi (or by using ^{qUi)) = 
fiiiqU^T) = liiU*q*)). 

We next show that (|12p and its equivalent for a = — 1 characterize 
uniquely ^ £ M when the tij are small enough. It will be crucial here 
that fi{q) is bounded independently of the tj's (here by the constant 1). 

Now, let iJ,,fi' £ M. be two solutions to SD[V,r] and set 

A{£)= sup \fi{q)-fi'{q)\ 
deg{q)<e 

where the supremum holds over monomials of C((f/j, [/*, ylj)i<i<m) with 
total degree in the Uj and Uj less than i. Namely, if the monomial (or word) 
q contains times Uj and aj times U* , we assume X^^i (o^ + ^7 ) — ^• 
Note that by traciality of /x, 

(13) Aii) = max sup |/i(C/f g) - /x'(C/f (/)| 

ae{+T,-i} degg<£-l 

and that by (I12p . we find that, for q with degree less than i — 1, 

HU*q)-^,'{U:q)\ < Y l(/^-^')(ft)l+ E l(/^-/^')(92)| 

q=qiUiq2 q=qiUiq2 

+ Y \if^-f^'){U:qi)\+ \if^-f^'mq2)\ 
q=qiU*q2 q=q\U*q2 

+ X]*iiK^ - ^J''){q^jq)\■ 

j 

A similar formula holds for \fi{Uiq) — fi'{Uiq)\ by conjugation, and therefore 

e-2 e-i 
A(^) < 2 ^ A{p) ^(P) + nDeA{e + D - 1) 

p=i p=i 

where we used that deg(gi) G {0, 2}, deg(g2) G {0, 2} (but 

A(0) = 0) and deg(gjj) < D and assumed \ti\ < e. Hence, we have proved 
that 
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£-1 

A(^) < A'^Aip) + nDeA{i + D). 
p=i 



^(7) < T^Hij) + ^H{j) 



Multiplying these inequalities by 7^ we get, since H{'-f) := Yle>i 7^A(i?) is 
finite for 7 < 1, 

7 , nDe 

1 — 7 7-^ 

resulting with -^(7) = for 7 so that 1 > + Such a 7 > exists 
when e is small enough. This proves the uniqueness. 

□ 

As a corollary, we characterize asymptotic freeness by a Schwinger-Dyson 
equation, a result which was already obtained in [32], Proposition 5.17. 

Corollary 3.2. A tracial state fi satisfies SD[0,r] if and only if , under fi, 
the algebra generated by {Ai, 1 < i < m) and (Ui, U* ,1 < i < m) are free 
and the Ui 's are two by two free and satisfy 

^([7») = Va G Z\{0}. 

Proof. 

By the previous theorem, it is enough to verify that the law fi of free 
variables {Ai, Ui, U*)i<i<m satisfies SD[0,t]. So take P = C/f/^i • • • C/f/Bp 
with some B^s in the algebra generated by {Ai,l < i < m). We wish to 
show that for all i G {1, • • • , m}, 

H ® nidiP) = 0. 

Note that by linearity, it is enough to prove this equality when iJ.{Bj) = 
for all j. Now, by definition, we have 

Ik 

diP= Yl E K'B^ ■ ■ ■ Bk-M ® Uf-^'Bu ■ ■ ■ u^;Bp 

k:ik=i,cik>0 1=1 



Taking the expectation on both sides, since fJ-{Uj) = and tJ-{Bj) = for 
all i 7^ and j, we see that freeness implies that the right hand side is null 
(recall here that in the definition of freeness, two consecutive elements have 
to be in free algebras but the first and the last element can be in the same 
algebra). Thus, iJ,{diP) = which proves the claim. 

□ 
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4. Formal solution and analyticity 

We have shown in Theorem 12.21 that the hmit points of the formal model 
also satisfy an equation similar to Schwinger-Dyson's equation. The only 
difference is that one of these equations is on the space of tracial states while 
the other one is on the space of tracial power states. In order to prove that 
the formal model matches the matrix model we need to study this formal 
equation and show that the series have a positive radius of convergence, 
hence providing a solution to SD[V,r] as defined in Definition 13. 1[ 

Definition 4.1. Let Vt = '^i'tili be a polynomial. Let t he a tracial 
power state in -A^|(Ai)i<i<„ • ^ tracial power state £ Ai is said to satisfy 
Schwinger-Dyson equation SD-^[T45'7"] if and only if for all P E C((j4j)i<j<m), 

m(P) = t(P) 

and for all P G C((f/i, U* , Ai)i<i<rn) , all i G {1, • • • ^m}, 

fi fi{diP) + fi{D,Vt P) = 0. 

(Here, both terms of the above equality are elements o/C[[t]] and the equality 
is formal.) 

We already know, due to Theorem 12.21 that there exists a solution to this 
equation. We now prove that this solution is unique. 

Theorem 4.1. There exists a unique tracial power state t ^ fit which 
satisfies Schwinger-Dyson equation SD-^[Vt5''"]- 

Proof. 

Let ;Ut be a tracial power state solution of SD-^[Vt,T]. There exists a fam- 
ily = {h,--- , kn) G N" in the algebraic dual of C{{Ui, C/*, ^i)i<i<m) 
such that for all P, 

keN"i=l *■ 

We will now show that the /^k are uniquely inductively defined by the relation 
given by SD-^[Vt?''"]. Let us define Ij the vector in which vanishes on 
every coordinate except the j-th which is 1. We get the following equalities, 
for all k, 

(1) If P is in C((^,)i<*<m), IJ'HP) = T(P)lk=o, 

(2) If P = RUiS with S in C{{Ai)i<i<m), f^^{P) = iiHSRUi), 

(3) If P = RU*S with R in C((Aj)i<j<m) and S does not contain any 
Uj (but may contain the ^7*), fi^iP) = fi^{U*SR), 

(4) If q does not contain any Uj, 

I^Hu:q) = - (t) E l^'''{U*q,)fi^"{U*q2) 

q=qiU*q2 ^ ^ k'+k"=k 

+ Yk,fi^-^^{UtqD,q,). 

i 
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(5) And for all q, 

g=gi(7,g2k'+k"=k ^ ^ 

+ E E (k/)/^'''('?i)^''"fe)-E%^''"'^(^'9^-^f^)- 

g=qiC/*g2k'+k"=k ^ ^ j 

One can see that this allows to compute uniquely any fi^{P). The first 
relation takes care of the non random case, the relations 2 and 3 use the 
traciality to place a variable [/ in a convenient place. Finally relations 4 and 
5 allow to compute lj}^{P) as a function which depends on the p)^ {Q) with 
degQ < degP and k' ^ k (first terms) or on the /z*^ {Q) with k' < k (last 
term). This is a well founded induction. Thus the y}^ are uniquely defined. 

□ 

We next show that this solution is not only formal but that it gives a 
family of solutions of the non- formal equation SD [Vt , r] , which depends 
analytically on the parameters (ti)i<j<n- 

Theorem 4.2. There exists e > such that for t S C",maxi<j<„ < e, 
the formal solution jj,t o/SD^[Vt,r] is indeed a convergent series. For all 
polynomials P, t £ B{0,e) = {t £ C" : maxi<j<„|tj| < e} — > fitiP) is 
analytic. 

In other words, there exists a family f/u'^,k = {ki, ■ ■ ■ ,kn) G in the 
algebraic dual of C{{Ui, U* , Aj)i<j<m) such that for all P, 

keN"j=i *' 
converges absolutely for maxi<j<„ \ ti\ < e. 

An immediate consequence of this result is to deduce that the formal 
solution is a real solution of SD[Vt,r] in a small parameters region, and 
therefore by Theorem 13. H equals the real solution. This will be a key to 
prove Theorem 10.11 (see section [7]) . 

Corollary 4.1. For small t, the formal solution of Schwinger-Dyson equa- 
tion SD^ [Vt , r] converges as a series. In addition, it matches the real solu- 
tion o/SD[Vt,T] which thus depends analytically in the parameters t of the 
potential in a neighborhood of the origin. 

Let us now prove Theorem 14. 2[ 
Proof. 

According to the proof of Theorem 14.11 the fi^ are uniquely defined by 
the family of relations (l)-(5). We only need to control the growth of the 
coefficients fJ'^iP) to show that Ht{P) is indeed convergent for small enough 
parameters. 
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To bound these quantities, we use the Catalan numbers 
Co = liCfc+i = ^ CpCk-p 

and the fact that they do not explode too fast; Ck+i ^ 4Cfc. We denote 
Ck := Hi Ck, and Dk := A''-'^Ck-i for k > 1, Dq := 0. The two key prop- 
erties of this sequence is first that it is sub-geometric (Dfc+i ^ iAD^) and 
secondly it satisfies = A Ylo<p<k ^v^k-p- Now our induction hypothesis 
is that there exists A,B>0 such that for all k, for all monomial P of degree 

(14) ^ C^B^Dp. 

We prove this bound by induction, and the relations (l)-(5) which define 
the fj}^. For k = (0, ••• ,0) this bound is satisfied since Dp > 1. We will 
check the induction for a polynomial of the form qUi since it is the most 
complicated case. 

IfJ'HQm ^ ^ If^""' {qiUi)\ 



k! ^ k'! k"! 

k'+k"=k 



^ k'! k"! ^ (k-1,)! 

'3=91 C/* 92 kj^O ^ ■>' 



k'+k"=k 

Now we use the induction hypothesis. If q is of degree p — 1, 
\l^HQUi)\ ^2 C^'B'''DgC^„B''"Dp_, 



klC^B^Dp - C^B^R 

k'+k"=k 



C^B^Dp 



T-r Ck.+i 1 (4A)^ 

i * 

The point is that we can choose A,B>0 such that this last quantity is 
lesser than 1. For example take A > 4"+^ and then B > 2nD{4:A)^ . 

Thus, for ||t|| := maxj < 1/4S, for ah P in C{{Ui,U* , Ai)i<i<m) , the 

series X^^rii tif^ (B) is absolutely convergent. 

□ 



5. Combinatorics. 

The purpose of this section is to provide a graphical approach to the solu- 
tion of the Schwinger-Dyson equation, and therefore to the computation of 
unitary matrix integrals and free entropy (see sections El [7| and [8]) . Actually, 
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the proof of Theorem 14.11 gives a recursive way of computing formal solu- 
tions to the Schwinger-Dyson equation, and therefore numerical solutions 
with arbitrary precision. 

Before giving a detailed description of our combinatorial model, we start 
with an overview. We need the notions of a star, which is a pictorial en- 
coding of a monomial of C((C/j, C/*, Aj)i<j<m), of root star, which is a 
distinguished star, and of a map, which is a specific planar decoration over 
a set of stars and one root star. 

The goal of this section is to show that the limits of integrals on the space 
of unitary matrices are generating function of the number of some maps 
as described above. However we are not interested in all maps, but rather 
on some that arise from an admissible construction, which leads us to the 
concept of admissible maps. Last, we need the notion of weight of a map, 
and our result will be in terms of sum over admissible maps of weights. 

Let us point out that for the sake of clarity, although our natural play- 
ground is the algebra C{{Ui,U*,Ai) 

i<i<m) and our definitions work in full 
generality, we restrict ourselves in the examples to the case of one single 
unitary matrix U and two variables Ai =: A and A2 =: B. We first start 
with the definition of stars and root stars, in the spirit of [I6t [T7]. 

Definition 5.1. (1) A star is a circle endowed with the clockwise ori- 
entation, decorated with elements such as colored incoming or out- 
going arrows, and colored diamonds. One of the element is marked. 

(2) To each letter X-i in the alphabet {Ai,Ui,U*)i<_i<m, we associate 
(bijectively) an element as follows; a diamond of color i if Xi = Ai 
and a ring of color i if Xi = Ui or U* ; in the case ofUi (resp. U* ) we 
attach before the ring an outgoing arrow of color i (resp. we attach 
after the ring an incoming arrow of color i ) outside of the circle. 

(3) To a monomial q £ C((C/j, f/*, ^i)i<i<m)j we associate in a canonical 
way a star of type q by drawing on the clockwise oriented circle the 
elements associated to the successive letters of q, while the element 
corresponding to the first letter of q is marked (or distinguished). 

(4) A root star of type q is obtained by drawing on the oriented circle 
the elements associated to the successive letters of q in the counter 
clockwise order, the arrows being drawn inside the circle. Its first 
element is distinguished. Although the maps are on the sphere, in 
the graphical representation of this section we will draw them on 
the plane and, to highlight the role of the root star we will draw it 
in this section such that it contains all the other stars. It can be 
viewed as the star centered in infinity or as the outer face of the dual 
map. Besides, on a root star we will distinguish a root element. If 
q contains no Ui nor U* , there are no root element. If q contains 
a Ui, the ring associated to the last {Ui, 1 < i < m) is the root 
element. If q contains no Ui but some U^ , the ring associated to the 
first ([/*, 1 <i <m) is called the root element. 
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(5) A multistar is a set of k stars inside a root star drawn on the same 
plane with a coherent orientation. 
The figure U\ shows a concrete example of a multistar. In the middle of 
the picture there is a star of type U*AUB and, surrounding it, a root star 
of type U*A^UB^U*A^UB. 



distinguished element 

root element 

A 
B 

oriented edge 
ring 

orientation 



Figure 1. Star of type U*AUB and root star of type U* A^U B'^U* A^U B . 

We are now ready to introduce the main objects in our combinatorial 
model, namely, maps: 

Definition 5.2. A map is a decoration of a multistar into a connected 
graph embedded in the plane by drawing two species of edges between rings: 

(1) A first category of edges, called "dotted edges", can be drawn between 
rings either attached to two outgoing arrows of the same color or to 
two incoming arrows of the same color. These edges can only have 
rings as end points, not diamonds or arrows. Rings can have any 
number of dotted edges going out of them, possibly none. 

(2) A second category of edges, called "colored oriented edge" arises from 
the connection of an arrow going out of a star (associated with a 
variable Ui ) into an incoming arrow ( associated to a variable U* ) of 
the same color. These colored oriented edges is a pairing between the 
set of Ui 's and the set of U* 's: exactly one incoming arrow is glued 
to each outgoing arrow. 

In addition, all the above edges do not cross, all arrows are paired but rings 
can be attached to any number of dotted edges (including to none). 

In the remainder of this section we keep considering pictures drawn on 
the sphere (and in fact on the plane). They therefore give rise to graphs 
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with vertices, edges and faces - together with additional decoration. For our 
forthcoming definitions, we need to clarify the notion of 'face': we consider 
that faces of a graph are the connected components of the complementary 
of the graph on the sphere. However, we take the convention that the 
original stars are 'fattened vertices'. Therefore the interior of stars will not 
be considered as faces (neither is the exterior of the root star). 

Each 'face' component of a map is isomorphic to a disc; thus this is an 
actual face. This is due to the fact that our map is embedded into a sphere. 
This condition would not be granted in the case of an embedding into a 
higher genus oriented 2D compact manifold. In this case it would have to 
stand in the definition of a map of 'higher genus': this will be of use for 
future work but for the sake of simplicity we do not emphasize this notion 
in this paper. 

Next, we define the weight of a map. The boundary of a face is homeo- 
morphic to a circle, it is given an orientation (the orientation of the sphere) 
and is decorated with diamonds (note that all arrows have been paired); it 
thus has the structure of a star except for the distinguished element. 

Definition 5.3. Assume we are given the tracial state t of 

• First we define the weight of the faces of a map. The boundary 
of a face have the structure of a star, i.e. it has the topology of 
a circle with some diamonds on it. We can therefore associated 
each of these boundaries with a monomial in the Ai 's, given up to 
cyclic permutation (or equivalently up to knowing its first letter). 
The weight of a face is the trace T{q) (which does not depend on 
cyclic permutations) of the monomial q associated with its boundary. 

• The weight of the map m, denoted by Mm{T), is the product of the 
weights of its faces times a sign given by —1 to the power the number 
of dotted edges. 

As we said before, not all maps will contribute and we need to define now 
the notion of admissible maps. Admissibility can be checked by an inductive 
procedure IP, which ressembles Tutte's surgery [29] and which amounts to 
check one after the other whether edges of the map are admissible. Once an 
edge has been checked, it is frozen and we continue by checking the other 
edges. 

Inductive Procedure IP : 

a- If the root star has no root element, then it can not be connected to 
any other star and hence the graph can not be a map unless there is no 
other star in which case the map is just the trivial graph with no edges. 

b- The root star has a root element which is associated to a Ui (resp. a 
f/*), for some i G {1, • • • ,m}. 

1-Then, we first check the admissibility of the dotted edges starting from 
this root element. We first consider the dotted edge which is the farthest 
from the arrow and declare it admissible if its other vertex is a ring of 
an outgoing (resp. ingoing) arrow and that there is no other dotted edge 
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attached to this ring which is farther (amongst the unfrozen dotted edges) 
from its arrow. Once this condition is verified, we freeze this dotted edge 
and the root element remains the root element. We check all dotted edges 
of the root element inductively. Once a dotted edge has been checked to 
be admissible, it is frozen and we go on checking the others. Once all 
the dotted edges of the root element have been checked, they are frozen 
and may separate the graph into subgraphs. Thus, the map may have 
been cut into disjoint subgraphs whose boundary (which may contain dotted 
edges) is homeomorphic to a disc (In the case where it has edges glued with 
an internal star, we see these other stars as part of the external star by 
following all the graph connected to the external boundary). In each of 
these subgraphs, we declare the first (following the orientation of the plane) 
element (corresponding to a f/j or a U*) after the last frozen dotted edge 
of its boundary as distinguished. We then define the root element of the 
boundaries of these subgraphs by the same procedure as for the root star. 
The boundary of each subgraph is then a star and these subgraphs have now 
the structure of a map; we will call them submaps. 

For instance, in figure [21 once the two dotted edges have been checked 
and frozen, the map is cut into two disjoint submaps, the left one having a 
fixed frozen dotted edge between the root element and the inner star (thus 
forbidding other edges to cross it and allowing us to consider the inner star 
as part of the external star) . The boundary of this left subgraph is now seen 
as a star of type q = U* A^U BU* AU B . This left subgraph has the same 
distinguished element as before but a new root element (here the outgoing 
arrow on its boundary corresponding to the first U in q). 

2- When all dotted edges are frozen, we check that the arrow of the root 
element is paired with an arrow of the opposite direction (note that if the 
root element comes from a U* , it can only be paired with an element of 
another star since by definition there is no more outgoing arrows on the 
root star). The oriented edge is seen as a fat edge. In particular, if the 
oriented edge link the root star with another star, we see this other star as 
part of the root star for the next step, i.e we identify the root star of type 
QUiP glued to the star of type RUIS (by the marked Uis) with the root 
star of type PQSR with, by convention, the distinguished element chosen 
to be the closest element after the glued U* . If the oriented edge link two 
rings of the root star, two disjoint subgraphs are formed and we proceed as 
in 

c- We continue the inductive procedure on the submaps until all edges 
have been checked to be admissible. 

Now we can define weighted sum of admissible maps. 

Definition 5.4. Assume we are given the tracial state r o/ ([2]). 

We define the weighted sum of admissible maps constructed above the 
stars ri, • • • , r„ and the root star P: 

Mn,-,r„(P)=5^M„(T) 
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where the sum runs over all admissible maps m constructed above ri, • • • , r„ 
with root star P. Assuming that Vt = t\q\ + . . . + inQn where qi are mono- 
mials, we define the formal series: 

keN" 

with Mfc^^... ^fc^(P) = Wiq-^^... ^q^^... ^q^^... ^q^{P) whcrc the monomial qj appears 
in kj successive position and = Yi^i ' =Y[^i- 

Remark that we do not count all the maps which contain the stars ri,. . . ,r„ 
but only those that are constructed using our inductive rules; they for in- 
stance forbid to glue the two same rings more than twice. 

However, a given map is counted at most once since there is only one way 
to decompose it using the procedure IP. Indeed, it is easy to check that 
at each step we have only one possibility for the next step since the dotted 
edges have to be drawn one after the other following the orientation and no 
new dotted edge can be drawn after the arrow of the root has been glued. 

Example 

Let us show some examples. We start from one root star and a star on 
the sphere (see figured]). We want to construct maps above these stars with 
our rules, starting with the root element shown by the arrow outside the 
root star. Figures O [3] and [5] are examples of such maps. Note that the 
weights of the maps of figures [2] and [3] are the same, the only difference is 
the way the three rings are glued. There is a third way to glue those three 
rings shown in figure H] which is a map but can not be obtained by our rule 
of construction (and thus is not admissible). 




Figure 2. A possible map. Its weight is t^^{A^ (g) _B (g) i?^ (g) 

We now come to the main theorem of this section, namely the graphical 
expansion result for Mt: 
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Figure 3. Another one. Its weight is t'^^{A^ ® B ® B"^ ® ® B) 




Figure 4. A counterexample: IP is violated because of the 
order of the dotted edges at the root element 




Theorem 5.1. Let V = X]i<i^n^«^« ^ polynomial. Let be a solution 
of SD [Vt , r] and Mt be the formal series defined for monomials P by 

" t^' 

Mt(^)= E n^Mk(p) 

k6N"i=l 



26 BENOIT COLLINS, ALICE GUIONNET, AND EDOUARD MAUREL-SEGALA 

where Mk(i-*) is the weighted sum of planar maps with one root star of type 
P and ki stars of type qi. If we extend the definition of Mt hy linearity 
to any polynomial P then the series Mt(P) is absolutely convergent in a 
neighborhood of the origin and, 

Mt(P) =/it(P). 

Proof. 

For the sake of clarity we first prove the case ^ = and show that 
M(P) := Mo(-P) = /Ut(P) for a monomial P. 

We proceed by induction on the total degree in Ui, 1 < i < m, in q. 

Suppose that there is no variable Ui in P. Then either there is no variable 
U* and both sides of the equality are equal to t{P), or there is a U* and 
both sides vanish: the left hand side by freeness between Ui and the Ai's 
and the fact that all non-trivial moments of U is and the right hand side 
because one can not glue the arrow coming out from this C/* anywhere. 

We assume our identification proved when the degree of P in the C/j's is 
less than k. We next take q with degree in the C/j's equal to + 1. Thus 
we can assume that there is a U in P, and we consider the last one in P so 
that P = pUib with b a polynomial in the U* and the Aj^s, 1 < j < m. By 
definition, M.(j)Uib) = M(6pC/j) since it depends only on the position of the 
last Ui. Thus, we may assume that P is of the form QU with Q of degree 
k. We apply Schwinger-Dyson equation to this quantity: 

(15) ^i{QUi) = - KRUi) ^ t^isUi) + Yl KR)(^t^iS) 

Q=RUiS Q=RU*S 

Now, we can apply our induction hypothesis since all polynomials appear- 
ing in the right hand side have degree strictly smaller than k + 1. 

We need to show that this is exactly the induction relation for maps. To 
construct a map above a star of type QU, we first look at the root element 
U and we have to decide what to do first with the dotted edges. There are 
two possibilities: 

(1) The first possibility is that there is no dotted edge going outside of 
the ring of the root. In such a case, we can glue the arrow to any 
other arrow of opposite direction and of the same color (correspond- 
ing to a variable U*). This implies that Q decomposes into RU*S 
and we construct an oriented edge between Ui and U* . Thus we 
separate the map into two parts and we have to construct a map 
above the R part and another one above the S part (this is the case 
2 of IP). This gives 

M(i?)M(5) 

possibilities which is exactly the possibilities counted by the second 
term in the right hand side of (llSp . 

(2) The second possibility is that we glue the root ring to another ring 
with a dotted edge. Thus Q must decompose into RUiS and the 
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creation of the dotted edge amounts to decompose the map into 
RUi and SUi and again to continue the construction of the map we 
will have to construct a map above the RUi part and another one 
above the SUi part (note here that when a dotted edge is attached 
to a circle of an Ui, the arrow and the circle keep their structure and 
live on the right of the dotted edge). In this procedure, we have fixed 
one dotted edge and thus multiplied the contribution of the resulting 
map by —1 (this is the case 1 of IP). The resulting contribution to M 
is therefore —M.{RUi)'M{SUi). Thus, the first term in (fT5]l computes 
the operation of gluing rings by dotted edges. 

Putting these two possibilities together we see that the state n and the 
enumeration of maps M satisfy the same induction so that they are equal; 
M.{pUib) = fj,{pUib) for any b monomial which does no contain any of the 
(Ui, 1 < i < m). Note here that no dotted edges between rings of incoming 
arrows can be drawn since if there are no outgoing arrows in a map, but some 
U^, there is no contribution. By traciality of /i, we deduce as well that Mq is 
tracial. Indeed, if we decompose p, q into p = PiUi^p2Ui2 ■ ■ ■ Pn-iUi„_^Pn and 
q = qiUj^q2Uj^ ■ ■ ■ qr-iUj^_-^qr with monomials Pi,qi which does no contain 
any of the {Ui, 1 < « < m), 

M{pq) = M{ipqiUj,q2Uj^ ■ ■ ■ qr-iUj,._^qr-iUj^_,qr) 

= KPQiUj^Q2Uj2 ■ ■ ■ qr-iUj^^^qr-iUj^_^qr) = ix{pq) 
= f^iQP) = f^{{qPlUiiP2Ui2 ■ ■ ■ Pn~l)Ui^_^Pn) = M{qp). 

Now we turn to the general V case. 

We first check the induction relation when the root star P contains a Ui 
for some i G {1, • • • ,m} so that we can write P = QUi. Let us denote for 
n-tuples k= {ki, - ■ ■ ,kn) and i = {h, - ■■ , Z„), (^) = J], (/')• check the 
formal equality by considering the induction relation, now given by: 

f.'^^'KQu.) = - E C"^/')^HRu^)^f^^^'^-Hsu.) 

£<k+lj Q=RU^S ^ ^ 

(16) + E E (^V')AR)^^''^'^-'is) 

^<k+lj Q=RU*S ^ ^ 

- kjfi^iQUiSRUi) - kjfi^iQSR) 

qj=RUiS qj=RU*S 

We need to show that the enumeration of maps satisfies the same relation. 
We start by putting stars of type (g^, 1 < j < n) inside a root star of type 
QUi and we wonder what happens to the root element Ui. We apply one step 
of IP. Two things can happen. Either we link Ui to another part of Q and in 
that case we have already shown that the possibilities are enumerated by the 
first two terms of the induction relation. Here, note that the product of (^j') 
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corresponds to the possible distribution of stars in each part (or submap) of 
the map, since aU the stars are labeled. 

Thus we need to show that the two other terms take into account the case 
where Ui is linked to another star of type qj. According to our construction 
rules we have two possibilities: 

(1) Starting from Ui we glue the arrow to an arrow of the same color 
entering a star of type q. This rule forbids any other gluing from Ui, 
this is counted by 

^ k,f,iQSR). 

gj=RU*S 

The coefficient kj counts the number of choices for the star of type 
qj since they are all labelled. 

(2) The other possibility is to glue the ring to a ring of the same color. 
This leads to 

- kjfx{QUiSRUi) 
gj=RU,S 

possibilities. 

In the case where P does not contain any Ui, 1 < i < m but still some U*, 
the root of the root star can only be glued by a dotted edge to any other 
U*, or by a directed edge to a Ui of a star. The resulting induction relation 
is exactly given by the formula obtained by conjugation of p6|) . hence again 
-^k(-P) = 1^{P)- This completes the proof. 

□ 

This theorem gives a combinatorial interpretation in term of maps to the 
unitary integrals. The fact that we do not take the sum on all maps but 
only on admissible ones makes this interpretation less transparent than the 
one for the gaussian case found in [7j. However, now that we know that the 
series can be identified to the matrix integral, we obtain some combinatorial 
identities which show that IP is less rigid than it looks like. 

Corollary 5.1. Let V = Yli'^ili- 

(1) For all P, Q, 

Mt(PQ) =Mt(QP). 

(2) For all monomials ri, . . . ,r„,r„_|_i, and all permutation a of n + I 
elements, 

(3) Assume that we define another procedure to define the root element 
of the root star (for example we pick the root element to be the second 
ring avaible if possible, or we pick a ring at random, or any other 
choice which may change during IP for the root stars that are created 
during the procedure when new faces are added). This will change 
the notion of admissible maps and we can define a new weighted sum 
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^'n,--- ,r„{P) ^'^^ ^ '^s^ series M[{P) where the sum occurs on these 
new maps. For all ri, . . . ,rn, P, 



m. 



-ri,--- ,rn 



m: 



Mt(P) 



m;(p). 



Note that due to the definition of admissible maps via the procedure IP, 
those properties are far from being obvious from a purely combinatorial point 
of view. Still they will appear as an easy consequence of the identification 
with the matrix model. 

Obviously different roots lead to a different procedure IP, and thus po- 
tentially to different maps. It is actually possible to see through examples 
that this phenomenon actually happens. 

However, it follows from the second point of the corollary that the choice 
of the root does not affect the weighted sum. The first and third points 
show that the choice of the root element and of the root star does not affect 
the final series. We were not able to give a more direct combinatorial proof 
of that result. 

To be more specific on the impact of the choice of the roots on the maps, 
let us call clusters the equivalence class of rings for the equivalence relation 
generated by a ~ 6 if the ring a is glued to the ring 6 by a dotted edge. 
Changing the choices of the roots will lead to different admissible maps 
since it will allow different positions for the dotted edges. For example, they 
were three choices for the starting root in figured) For each of these choices, 
two of the three maps represented in figures [21 [3] and U] would have been 
reachable by the inductive construction IP but not the third one. The one 
who is not constructible depends on the choice of the first root. It seems 
that if the maps are different, nevertheless the clusters are the same and in 
that simple case, knowing this cluster is sufficient to define the faces created 
by the dotted edges and thus the weight of the maps. 



Changing the root element of a star is the same thing than making a cir- 
cular permutation of the variable of the associated monomial. The theorem 
shows that weighted sums are equal to the limit of the empirical measure of 
the matrix model which are tracial. The first and third items are a direct 
consequence of this identification. 

For the second item, observe that permuting the first n monomials doesn't 
change the sum by its definition. Thus we only need to show that 



Proof. 



M. 



-ri,--- ,r„ 



(P)=Mp,,,,..,,Jri). 



Let us define V = Uiri + tP. We will again use the identification with 
the matrix model but now we will use the formal version. The coefficient 
M^i,... ,r„(-P) appears as the coefficient of the limit tracial power state //-^ by 
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Corollary 13.11 and Theorem 15. II More precisely, 



QU 



/(P) 



Ui=0 



N dui ■ ■ ■ dUr. 

We now use the fact that /i-^ is the limit coefficientwise of the formal model 
defined in Thus, 



u,=Ot=0 



M,;=0,t = 



TV dtYl^duiN^ 

We conclude by noticing that this last expression is symmetric in the mono- 
mials ri,. . . ,r„, P. 

□ 



6. Application to free probability 

In this section we look at applications of the combinatorial results of 
section [5] to free probability. 

Let us assume that the C/j's are chosen independently according to the 
Haar measure. If we define Xi = U* AiUi then the XiS are asymptotically 
free (according to a theorem of Voiculescu [31]) and with fixed distribution 
/i uniquely defined by the distribution of the Ais. We are interested in using 
our setup to compute limits of moments of these variables or in other word 
to compute the moments of free variables: 

^{Xi^...Xi^). 

According to our interpretation this can be computed by looking at the maps 
above the star of type Xi^...Xi^ without any other stars, in other words we 
have to focus on computations of M((7) = Mo(g) which turns out to be equal 
to ^(g) where is the free state product (see Corollary 13. 2p 

We are interested in using this method to compute some non-commutative 
moments of free variables, in relation with Speicher's non-crossing cumulants 
theory, cf [27] . 

6.1. One star maps. For these purposes we need to find a simplified in- 
terpretation of M(g) in the single star map. 

For this case with only one star, the combinatorial interpretation can be 
slightly modified. First, we do not need to consider dotted edges between 
incoming arrows since if there is a there must be a Ui which can be chosen 
as the root element or we can not build any map. But the main difference 
is that now each time we glue two rings, the edge newly created separate 
these two rings into two different faces so that they can no longer be glued 
together. Thus, we can forget about the restriction of the construction rules 
and present a simpler description in that case. Instead of gluing the ring two 
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by two we will now glue them together. We define a new structure which we 
will call a node and now rings can only be glued to node and a node can be 
glued to any number of rings. A one star map is a map with one star where 
the arrows has been glued two by two while respecting the orientation and 
rings may be glued to exactly one node, each node is glued to an arbitrary 
number of rings but at least one. Figure [6] shows the new representation 
of a one star map. The trick to go from the previous interpretation to this 
one is to glue together to a node all the rings that are in the same class 
of the equivalence relation generated by being glued. In order to compute 
the weight of such a map, observe that several maps give the same one-star 
map, but the weight is easy to compute since as we will see we only need to 
add a factor C^-i for each node of degree d. 



Figure 6. Reduction of a map on one star to a one-star map. 

Definition 6.1. A one-star map is a connected graph embedded on a sphere 
above one star and with some edges such that 

(1) Edges are drawn only between rings and must not intersect. 

(2) Arrows must be glued two by two while respecting the orientation and 
the color: an arrow going out of a star (associated with a variable 
Ui) is always glued to exactly one other arrow going into a star (as- 
sociated to a variable U* ) of the same color. This pair of arrows 
creates an oriented edge. 

(3) Any number of rings may be glued together on a node. 

The weight of a one-star map is the product of the weight of its faces 
which is defined as before as trace of products of Ai^s times the product of 
the weight of the nodes. The weight of a node of degree d is {—l)'^~~^Cd-i. 
We define Mo((7) the weighted sum of one-star map above a star of type q. 
Note that we no longer need to take care of roots and of maps that can be 
built with some set of rules. 

Proposition 6.1. For all monomial q, 




tx{q) = Mo(g). 
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Proof. _ 

We only need to show that M(q) = Mo(q'). For this we need to compute 
the number of maps above one star that are reduced to a given one-star 
map. The reduction goes as follows: two rings are glued to the same node if 
they are linked by a sequence of dotted-edges. We have to count how many 
configurations of dotted edges lead to a node of degree d. When one of the 
ring glued to this node becomes the root in the recursive construction, it 
has to be glued to one of the other ring glued to the node. Thus it separates 
the set of ring into two subsets, so according to our inductive procedure 
of section [5l we have to continue to glue this ring to other ones while we 
continue the construction in the face newly created. This yields a structure 
of tree on this set of rings. We have as many choices as they are trees with 
d — 1 edges (to glue the d ring we need exactly d — 1 edges). This explains 
the factor C^-i- The factor (—1)'^"^ simply comes from the factor —1 which 
comes with each edge. 

□ 

6.2. Maps and cumulants. Let Ai, . . . , An be self-adjoint variables and 
U a unitary matrix, free from the Ai^s. Then choosing k indices ii, . . . ,ik in 
{1, n} one has 

^i{Ai, ... A, J = fi{U*UAi, . . . U*UAi^) 

Let us apply Schwinger-Dyson equation with respect to U to the above 
equality, and let us rearrange the sum according to the non-crossing partition 
of Ai^s generated by the oriented edges. Obviously one obtains a formula of 
type 

(17) fi{A, . . . A,^) = K^(.Ai„...,A,^) 

■K&NC{k) 

where NC{k) is the non-crossing partitions and Kj^ is a /c-linear form mul- 
tiplicative along the blocks of vr in the sense of Speicher: if vr = {Vi, . . . , Vn} 
with the block Vi = {a\, . . . , a* .} 

KTr{Xl . . . Xk) = -^(r.,) (A'„» , . . . , Xa^^^ ) 
i 

where (rj) represents the partition on rj elements with only one block. 

The fact that such a formula holds true for any choice of non-commutative 
laws for ^j's proves via the moment-cumulant formula that K-,^ has to be 
Speicher 's non-crossing cumulants K-,^. But it is also given as a sum on maps 
by our graphical model. 

Let us recap this in the following proposition: 

Proposition 6.2. The n-th non-crossing cumulant of the variables Ai, . . . , 
Ap is the weight of all one-star maps over the star build by putting in the 
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clockwise order a ring, a diamond of color ii, a ring, a diamond of color 
i2, - ■ ■ , a ring, a diamond of color ip. 

Note that we have defined this map above a star which is not of type q 
for any monomial q. This would be a problem for admissible maps since IP 
requires the presence of oriented edges. But the definition of one-star map 
is fine in this context. 

Actually, Proposition 16.11 gives us a new proof of the following Corollary, 
due to Speicher and known as non-crossing Moebius formula 

Corollary 6.1. The following inversion formula holds true: 

neNC'in) B block of-n'^ 

where tt'^ is the Kreweras complement (see and Cq the Catalan number. 
Proof. 

This is a direct consequence of the previous proposition. Remember that 
Kn{Ai, . . . , A„) is a weighted sum over maps with dotted edges since the 
star contains some rings and no arrows. These dotted edges form a non- 
crossing partition of by saying that two rings are in the same 
component if their are linked to a same node. The weight associated to this 
map is a product whose factors are: {—l)'^~^Cd-i for each node of degree d 
and the weight of each face. The faces are by definition the component of 
the Kreweras complement of vr'. Thus we obtain: 

■K'&NC(n) B block of it' 

The formula follows after taking vr' = vr'^. 

□ 

As a further remark, one can also read graphically the main properties of 
cumulants, for example, Kn{Xi, . . . , Xn) = as soon as there are occurence 
of free elements. More precisely, assume that we can partition the XiS into 
two families the Aj^s and the B^s with the algebra generated by the ^j's 
free from the algebra generated by the -Bfc's. Then if all the XiS do not take 
value in the same algebra, Kn{Xi, . . . , Xn) = 0. Indeed, one can replace all 
the family of Aj^s by the one of V* AjV with V unitary and free from the 
other variables. Now when looking at the combinatorial interpretation of 
^{Xi, . . . we can see that the oriented edges coming from V separate 
the components containing the Aj's from the others. By following those 
edges we see that the faces they are defining contain only variable from 
one of the two algebras (The edges are going in the clockwise order around 
the faces which contain the B^s and in the counter-clocwise order around 
the faces containing the Aj's). Thus, in the decomposition (jlTh . the terms 
corresponding to partitions with one component containing both some AiS 
and some Bj^s vanish. By uniqueness of the decomposition into cumulants 
we deduce that those elements vanish i.e. Kn{Xi, . . . , Xn) = 0. 
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These remarks are not new but this shows that our graphical model fully 
encompasses the theory of non-crossing cumulants and that the Schwinger- 
Dyson equation can also be read in terms of cumulants. 

It is interesting to mention here that papers [24j and [25] have developed 
a calculus on annuli which seems to be related to our graphical model. 
However these approaches only deal with the asymptotics of second order 
cumulants whereas our approach via formal calculus, see section HI allows 
us to deal with arbitrary order cumulants. 

The actual relation can be found in [10] , where convolution on partitioned 
permutations is introduced and showed to be the relevant algebraic tool to 
handle higher order freeness, namely, the asymptotic behaviour of cumulants 
of unitarily invariant random matrices. 

But the results in our paper give an explicit algorithmic description of 
the Moebius inversion formula and therefore of higher order cumulants. As 
in the one star case, cumulants are also obtained by inserting an outer U*U 
between each variable of each star and by looking at generating function 
where U is linked to its neighboring U* . 

It is interesting to see that a direct (yet difficult to describe) graphical 
reading of the Schwinger-Dyson equation, which is our main tool of inves- 
tigation of unitarily invariant matrix models, yields non-crossing and could 
yield higher order moments related series and operations similar to convo- 
lution, although these latter results rely on more representation theoretic 
grounds (Weingarten function theory as developed in jllj). 

It is not obvious to us how the Schwinger-Dyson equation can be read 
off from the results of [lOj (without writing a change of variable invariance 
formula), and it would be interesting to attempt to figure out the meaning 
of Schwinger-Dyson equation at the representation theoretic level. 



7. Application to the asymptotics of lN{V,Af) 

In this section, we investigate the free energy by using the combinatorial 
interpretation of the previous section. 

Let {qi,--- ,qn) be fixed monomials in C((C/i, C/*, ^j)i<i<m), let V = 
YlitiQi be a self-adjoint polynomial and iNiViAi) be given by 

Theorem 7.1. There exists e = e(5i, • • • , q-n) so that for any t € C"nS(0, e) 
such that V = V* for any a S [— 1, 1], 

Fy,,(a):=Jim^^log/^(ayt,Af)= Yl fl^-^^^Ulir ■ ■ ,qn,T). 

keN"\(0,.. ,0)4=1 *■ 

Moreover, 

Mkiqi,--- ,qn,T) = Yl ^™(^) 

madmissible maps with ki stars qi 
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is the weighted sum of maps constructed above ki stars of type qi for all i, 
after choosing one of them as a root star (this is well defined according to 
Corollary \5.1\) . 



Proof. 

Let 



Then, if a G M, 

Assume that t is small enough so that Corollary 13.11 holds and remark that 
Vat is self-adjoint and such that \ati\ < e for all i and all < a < 1. Thus, 
for a e [0, 1], 

lim daF^^ = fiatiVt) 

with ^at the solution to SD[aVt)''"]- By dominated convergence theorem 
(since daF^t uniformly bounded in N and a € [0, 1]), we deduce that 



where we used that = 0. 



lim = / fJ.at{Vt)da 



□ 



Here also, we obtain the following important corollary, as a consequence 
of Corollary 14.11 



N^oodz'^ Jym 



Corollary 7.1. The following holds true: 

Qk 

Af— >oo _ 

Qk 

= Q^Fv,r{z)U=0 

In particular, this result allows us to give an expansion of the Harish- 
Chandra-Itzykson-Zuber integral as a generating function of the number of 
some maps. Let us recall the exact expression of this integral: 

F^'^'iz) := ^log HCIZizA,B) = ^log e^^^^^^'^^^Uu. 

The maps appearing in the expansion contain only stars of type U* AUB 
(see the star in the middle of figure [T|) . Besides we can build these maps 
without considering the rings attached to variable U* since we will always 
be able to choose the root element to be a f/ (a U* always comes with a U 
for this potential). 

Since the number of diagrams is growing quickly we compute only the 
first term of the expansion. Note that when gluing the arrow of the root of 
the root star, we must always glue it to another incoming arrow of another 
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star and hence we shall never see the case of a root star with no C/j's. Again, 
we therefore do not see dotted edges between incoming arrows. 

Besides, we consider only the case where the distribution is centered, that 
is when t{A) = t{B) = 0. The other cases can be deduced easily from this 
one since we have the relation 

pa+A,h+B^^^ = Fp^'iz) + ^{bTiA + aTiB) + zab. 

In terms of diagrams, this means that we only need to consider diagrams 
such that no face contains only one diamond. 

According to the previous theorem, limAr^oo Fn' (z) has, for small z, an 
expansion Yln^n^"'- We now use this graphical representation to compute 
the first terms of this integral. 

Since the distributions are centered, the first term Fi is zero. 

The second term F2 consists of maps constructed with two stars of type 
U*AUB. There is only one way to add edges between these two stars to 
construct a connected map without faces which contains only one diamond, 
this is represented by figure [71 We obtain a map with two faces. One has two 
diamonds associated to A and the other one two diamonds associated to B. 
Thus the weight of this map is t{A'^)t{B'^). Since there is no gluing between 
the rings they are no other signs. They are only one way to distribute the 
labels on this picture (that is the second distribution leads to the same map) 
thus to obtain F2 we only need to divide by 2!, 

F2 = ^t{aMb'). 




Figure 7. Second term in the expansion of the HCIZ integral. 

We can continue this for the next terms in the expansion, the third term 
(see figure [8|) is in the same spirit and leads to 

F' = ^riA^MB^). 

The fourth term is the first one where gluings between the rings appear. 
Thus weigths with negative coefficients can occur. The sign of a map is easy 




to compute, it is —1 to the power the number of dotted hues in the map. 
Equivalently since in the case of HCIZ integral the number of oriented edges 
is equal to the number of stars, this number is also equal to the number of 
faces of the map and thus to the number of factor in the product of moments 
of the weight. In figure [U we have drawn all unlabelled planar maps one can 
construct with 4 stars. To compute the exact coefficient of each map one 
has to multiply it by the number of way to distribute the labels and divide 
by 4!. 

This leads to, 

= \r{A')r{B') - \r[A'fr{B') - \r{A')r{B'f 

+]^r{A^fr{B^f + -^T{A^)\{B^f . 

Here the weight are given in the same order than the maps in the figure. 
Note a new and interesting feature that appears in the third map: two rings 
are linked by more than one dotted edge. 

The other terms can be computed in the same way, for example figure [TOl 
represents the fifth term and gives 

Fs = -t{A^)t{B^) - t{A^)t{A^)t{B^) - t{A^)t{B^)t{B'^) 
5 

Thus the first terms agree with the expansion given in [36] on page 23, 
besides this allows us to answer a question raised in this paper. Indeed, 
the authors ask if there is an explanation to the fact that the coefficient 
of Fn all seem to be integer multiple of -. This is easy to prove with this 
graphical interpretation. To compute the contribution of a given unlabelled 
map we must distribute the labels {1, .., n} on its stars, count the number of 
different map that we obtain and divide by n!. But after choosing the star 
which received the label 1 we have (n — 1)! ways to distribute the remaining 
labels and they all lead to different maps (note that on the other hand, due 
to possible symmetry in the unlabelled map, different choices for the star 
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Figure 9. Fourth term in the expansion of the IZ integral. 

with the label 1 may lead to the same maps). Thus the coefficient in front 
of this map is a multiple of ^^'l^' = 1/n. More precisely it is 1/n times 
the number of choices of the star which carry the label 1 that will lead to 
different maps, in particular it is always less than 1. 

To finish, we wish to point out that we can recover results in [8] and 
|15j about scalings of HCIZ integral. In these two papers, one considers the 
scaling where A has small rank, which amounts to considering only terms 
t{A^) X P{B). Here the transformation depicted in section [6] applies and 
one sees that P{B) has to be Kk{B). In particular this means in the case 
that ^ is a rank 1 projection, that A'^"^ log HCIZ tends to the primitive of 
Voiculescu's i?-transform. 

8. Application to Voiculescu free entropy 

Voiculescu's microstates free entropy is given as the asymptotic the vol- 
ume of matrices whose empirical distribution approximates sufficiently well 
a given tracial state. Up to a Gaussian factor, it is given by 

= limsup limsup-^log/i|™(rij(^,e,A;)) 
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Figure 10. Fifth term in the expansion of the IZ integral. 

with UN the Gaussian measure on Tiiy and Tr{ij,, e, k) the microstates 

Tiiifi,e,k) = {Xu--- ,X^enN. \^Ti{Xi^ ■ ■ ■ X,^) - fi{X,^ ■ ■ ■ XiJ\ <e 

p < k,ii e {!,■ ■ ■ ,1^}, \\Xi\\oo < R}- 
When m = 1, it is well known [30j that fi £ ■p(]R) and 

Xin) = = j j \og\x- y\dfi{x)dfi{y) - \ j x'^dfi{x) + const. 

Moreover, one can replace the limsup by a liminf in the definition of x- 
Such answers (convergence and formula for x) are still open in general when 
m > 2 (see [5] for bounds). However, if is the law of m free variables with 
respective laws /Xj, then these questions are settled and 

m 
i=l 

We here want to emphasize that our result provides a small step towards 
dependent variables by showing convergence and giving a formula for the 
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type of laws fi solutions of Schwinger-Dyson's equations SD[V,r]. Indeed, 
we shall prove that 

Theorem 8.1. Let fi be the law of m self-adjoint variables Xi with mar- 
ginal distribution (fii, • • • , ^m)- Assume that Xi can be decomposed as Xi = 
UiDiU* with Ui unitary matrices in such a way that the joint law v of 
{Di,Ui,Ui)i<i<m satisfy SD[V,r] with r the law of m free variables with 
marginal distribution fj-i,--- , o^nd some potential V = Y17=i^i1i- 
sume that the ti's are small enough so that Corollaru \ 3.1\ holds. Assume 
also that the hypotheses of Theorem \7.1\ hold. Then, 



Xip) = lim inf liminf log/i|™ {TR{^l, e,k)) 

and a formula of xip) can be given in terms of the fi^ 's of Theorem \4.^ 
Proof. 

Indeed, let us consider V = V{UiAiU* , 1 < i < m) with V a self-adjoint 
polynomial and fi the unique solution of SD[V,t] with r the law of the 
Ai,l < i < m which is now chosen to be the law of m free variables with 
marginals distribution fii, 1 < i < m. Under the law /i^™", we can diagonal- 
ize the matrices Xi = UiDiU* with Ui following the Haar measure on IAm, 
and id d is the Dudley metric, we find that for sufficiently large 

= [d{fiD^, Hi) < e;Ac/iAC/*,i<j<m e rij(/x, e. A:)) 

= .,<A I du^---duA n ^^^(^*) 



where we denoted A(Aj) = Ilfc^j — Xj\ and daj\f the probability measure 
dar,{X):=Z-/ll\X,-X,\'e-^^^''^' J] 

In these notations, Di = diag(A|, • • • , A^) and A = (Ai, • • • , Xn)- Hereafter, 
uj^p i denotes the empirical ditribution of {£'i}i<j<„; u^p i (P) = 

N^^Ti:{P{Ei, 1 < i < n)). As a consequence, applying the large deviations 
result of j3j to the diagonal matrices Di, we find that there exists o(l) going 
to zero with e such that 

L^<e^'&^(A'0+^'°a) sup / dU^---dU^ 

l|Oil|oo<-B 
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with for k greater than the degree of V, 

L]v= sup / e^Tr(y)-^Tr(v)^^^...^f^^ 



l|Dil|oo<-R 

<*(Ag.,M,)<e 
l|ojoo<fl 



Now, for fixed R, any Di,D'- in d{fi^,,fii) < e, ||-Dj||oo < 



^ log In{V, Di) - ^ fog In{V, D[] 



with ?7(e, R) going to zero as e goes to zero for any fixed R. Hence, 

nmsup-^log/iv(V^, A) < F{V,fii) +r]{e,R) 

with F{V, fii) the hmit of N'^ fog/Ar(V", Ai) given in Theorem 17. II when the 
distribution of the Ai converges to free variables with marginal distribution 
fj-i. We thus have proved, letting e going to zero and then R,k to infinity, 
that 



X 



ifi)<f2lM-KV) + F{V, 



4 = 1 



Conversely, we have 



L 



AT 



N 



with 

Li := inf / dUi ■ ■ ■ dU^, 

-■~N 



I 

^ y^i^i^l )\<i<m 



1 • • • "<^m 

\\Di\\oo<R 
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for any 6 < e. Now, choosing 6 and using the continuity of u^rr ^ rr.i in 

\UiUiU ^ j'l<i<m 

the distribution of the uniformly bounded variables Di, we find by Corollary 
13.11 and our hypothesis that 

lim mr = 1 

Je^^''(^)dUi---dUm 

which insures that 

m 

1=1 

Thus we have proved that 

m 

x{^l) = Y.I{^li)-^l{v) + F{v,^l^). 

i=l 

Note that ^liV) and F(y,fii) can be written in terms of the fi^ of Theorem 
2] by Theorem [m 

□ 



9. Generalization to integrals over the orthogonal group 

In a recent article [35j , Zuber shows that the large asymptotics of two 
matrix integrals (the integral with external magnetic field and the Harish- 
Chandra-Itzykson-Zuber integral) enjoy a universality property in the sense 
that they are the same (up to a proper rescaling) if one integrates over the 
unitary or the orthogonal group. This property was also obtained (but not 
explicitly stated) in the case of the Harish-Chandra-Itzykson-Zuber integral 
in [18] where the rate functions for the large deviation principle for the law of 
the spectral measure process of the Hermitian and the symmetric Brownian 
motion were shown to differ only by a factor two. The Harish-Chandra- 
Itzykson-Zuber integral is rather special in the family of angular integrals 
and one can compute many interesting related quantities, regardless of the 
group on which integration is taken (see [U |T3] ) . 

In this section, we generalize this universality property by relating the 
large A'^ limit of any small parameter integrals over the orthogonal group 
with its complex analogue. 

Let us define 

(18) lhiV,Af):= [ e^Tr(y(o„o*A^,i<.<-))^o^...^Q^ 

where {Af,l < i < m) are N x N deterministic symmetric uniformly 
bounded matrices, dO denotes the Haar measure on the orthogonal group 
Oj\f (normalized so that J^^^ dO = 1). In this section we will assume that 

y is a non-commutative polynomial in the Oj, 0*,Af^ with real coefficients. 
Here, O* = O* is the standard involution O*- = Oji. Observe that if P is 
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a polynomial, P{Oi,0*,Af,l < i < m)* = P* {Oi,0* , Af ,1 < i < m) so 
that we keep also the notation P* . 

We then claim that we have the following analogue of Theorem EH which 
shows that the first order of integrals over the orthogonal group is the same 
as on the unitary group (up to proper renormalizations) ; 

Theorem 9.1. There exists e = e{qi, • • • ,qn) so that for any t G M"nS(0, e) 
such that V = V* = tiqi, if we define 

1 



then Fy^ exists and 



I yr^ T-r tf' 



keN"\(0,..,0) l<i<n 

Moreover, 



Mkiqi,--- ,qn,r)= J2 ^'»(^) 

madmissihle maps with ki stars qi 

is the weighted sum of maps constructed above ki stars of type qi for all i, 
after choosing one of them as a root star . 

The proof is based on the fact that if ^y'^ denotes the law on O"^ given 

by 

M^^^m, • • • , dO^) := ^^^ef Tr(V(O.A^Af ,i<.<H)dO, . . . dO^ 

and fl^ is the empirical distribution of (Oj, O*, Aj, 1 < ? < m), then we have 
the analogue of Corollary 13.11 

Theorem 9.2. Assume that V = 'Yl, self-adjoint. Let D an integer 

and T a tracial state in -^|(yii)i<,<„ be given. There exists e = e{D,m) > 
such that if \ti\ < e, (i^ converges almost surely under fi^'^ to the unique 

2 

solution fit of the Schwinger-Dyson equation SD[V,r]. Moreover, p,j^' = 

2 

/ii ' (/i^) converges as well to this solution as N goes to infinity. 

2 

In fact, since then we know that fit{P) expands as a generating function of 
the Mk(gi, • • • ,gn;T)'s, Theorem 19. II follows readily since for any a G [0, 1], 

d^l-^\ogll,qVt,Af) = \ji^^{V) 

converges towards ^iit{V). 

Proof of Theorem 19.21 The proof follows the same lines as the proof of 
Theorem l2.lt we make the change of variables O = (Oi, • • • , Om) G 
^(O) = (^i(O), • • • , ^^(O)) G with 
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where the Pj are antisymmetric polynomials (i.e. P* = —Pj)- The only 
change is that now Pj{0) are matrices with real coefficients and the differ- 
entials hold in the direction of A]^ which are the antisymmetric matrices 
with real coefficients. For large enough, ^' is a diffeomorphism; it is as 
in the complex case a local diffeomorphism which is injective. As such, its 
image is open and compact. O"^ is not connected but the union of copies 
of SO%N) = {0 £ C'jv;det(0) = +e}, e = +1 or -1. Since det(^j(0)) = 

det(0j)det(e^^-'(O)) = det(Oj), ^ maps SO^^{N)xSO^''{N)x- ■ ■ SO^-^iN) 
into itself for each choice of G {1, —1}. Therefore, by connectedness of this 
set, ^'(S'O^i(A^) X • • • X SO'^"^{N)) is open and closed and therefore equals 
SO^'^iN) X 50=2 (AT) X . . . 50='"(A^). Thus, ^' is a diffeomorphism of O]}. 
Like in the proof of Lemma 12.11 we need to compute the Jacobian of this 
change of variable. The same arguments apply to show that 

I det J*(0)| = exp(^Tr$ + 0(1)) 

with $ the linear operator defined on antisymmetric matrices by 

^.A = Y,diPM- 

i 

A basis of A\j is given, for A; < by 

E\kl)rj = ^r=k,j=l -lr=lj=k _ 

Therefore, the trace of any linear endomorphism ip on A]^ defined by f{X) = 
^£ AiXBg, for uniformly bounded matrices Ai, Bi, is now given by 



Tr(v.) = Y.^iE\kir^{E\kl)))=^-^{Y^A'^,Bi,-Y,<Bi) 

k<l I k^l k^l 

= ^Y.^v{A'^Tr{B')+Tv{AeBj) 

I 

= -Y^Tv{A^)Tt{B^) + N0{1) 

since the operator norm of Ai and Bi is uniformly bounded, 0(1) is uni- 
formly bounded in A^. 

We can apply this bound to our case where Ai and Bi are given by 
diPi =: ^£ Ai (g) B£. The Ag and Bis are uniformly bounded since the Oj's 
and the Aj^s are and non zero for a finite number of £'s, thus we deduce that 

. m 

I det J*(0)| = exp(— TV ® TvidiPi) + 0(1)) 



with 0(1) bounded uniformly in N. Since 0(1) is uniformly bounded, we 
can now proceed exactly as in the proof of Theorem 12.11 to show that for 
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any r G {1, • • • , m}, 

As a consequence, for any limit point r of fi^ , any antisymmetric polynomial 

(19) r (g) T{drP) + T{DrVP) = 0. 

If P is symmetric, we claim that for any r G {1, • • • , m}, 

(20) r T{drP) = T{DrVP) = 

so that (fT9]) still holds. Indeed, if Q is a word in the (Oj, O*, Aj, 1 < i < m), 

Q = QlOrQ2 Q = QlO*Q2 

Q*=QiOrQ2 Q*=QiO*Q2 

= Yl QiOr®Q2- Y Qi®o;q2 

Q=Q*OfQl Q=Q*OrQl 

= Y io;Q2r^Qi- Y Qi^iQiOr)*. 

Q=QiO*Q2 Q=QiOrQ2 

Since the trace is invariant under transposition, we deduce that for all P, 
fi^{P*) = fi^{P) and thus, 

fi^ ^ fi^ {drQ + drQ*) =0. 

With the same method, we can deal with the cyclic derivative term. In- 
deed, since Dr{Q*) = — {DrQ)* , if we write V = Q + Q* , we obtain: 

fl^{DrV{P + P*)) = fl^{Dr{Q + Q*){P + P*)) 

= fi''{DrQ{P + P*)) - A^((A.Q)*(P + P*)) 

= fi^{DrQ{P + P*))- fi^{{P + P*)DrQ) = 0. 

To sum up, 

fl^{drP) = fl^{DrVP) = 

from which we get (|20p by going to the limit. Since any polynomial P 
can be decomposed as the sum of a symmetric polynomial {P + P* /2) and 
an antisymmetric polynomial {P — P*/2), we conclude by linearity that 
(jl9p holds for any polynomial P. By uniqueness of the solutions to this 
equation for sufficiently small parameters ti proved in Theorem I3.H the 
proof is complete. □ 
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